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Abstract 

We consider the titt scattering amplitude obtained by Colangelo, Gasser 
and Leutwyler (CGL), using dispersion relations and chiral perturbation theory to 
two loops. We point out that the input used by CGL for energies above 1.42 GeV 
is incompatible with a number of experimental results. Moreover, the error they 
adscribe to the 5^ phase at 0.8 GeV, an important part of their input, is much 
underestimated, and its central value is probably displaced. Then, we compare 
CGL results with direct fits to experimental -kit data; the outcome of this test is 
mismatch of the CGL amplitude, for several quantities by as much as 4 standard 
deviations. We also evaluate forward dispersion relations with the CGL ampli- 
tudes. We find that these dispersion relations are fulfilled less and less well as the 
energy increases towards 0.8 GeV. Moreover, the size that the experimental results 
for 7T7T scattering indicates for the two loop corrections in a chiral perturbation 
theory analysis is less than those required by CGL, which suggests that at least 
some of the CGL corrections are so large to cover biases in their nn amplitude 
and in the evaluation of the scalar form factor of the pion. We conclude on the 
necessity of repeating the analysis of CGL with correct high energy amplitudes 
and a more realistic input for the 5^ phase at 0.8 GeV. 



* This note contains the material presented in talks given at "Photon 2003" (Frascati, April 2003), "QCD at 
Work" (Conversano, June 2003), "Chiral Dymamics" (Bonn, September 2003), Benasque (July, 2004) and "Quark 
Confinement and the Hadron Spectrum" (Sardinia, September 2004). 
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-COMMENTS ON SOME CHIRAL-DISPERSIVE CALCULATIONS OF 7T7T SCATTERING- 



1. Introduction 



In two recent papers, Ananthanarayan, Colangelo, Gasser and LeutwylerW and Colangelo, Gasser and 
Leutwyler^ (to be referred to as, respectively, ACGL and CGL) have used experimental information, analyt- 
icity and unitarity (in the form of the Roy equations) and, in CGL, chiral perturbation theory, to construct 
the 7T7T scattering amplitude at low energy. 

We will mainly discuss the second paper, CGL. Here it is given a set of low energy parameters 
(scattering lengths and effective ranges) for S, P, D and F waves in which an outstanding precision (at the 
percent level) is claimed. For the SO, S2 and P waves 1 these authors also give a parametrization for the 
phase shifts valid for energies up to s 1 / 2 = 0.8 GeV. 

While the results of CGL constitute a substantial improvement over previous ones, the analysis of 
these authors presents a number of drawbacks. First of all, the input scattering amplitude which ACGL, 
CGL use for high energy (s 1 / 2 > 1.42 GeV) is not physically acceptable, as will be shown below; and the 
errors these authors take for some of their experimental input data are excessively optimistic. Finally, the 
input D2 wave used is incompatible with a number of requirements. Because of this, the CGL threshold 
parameters, and the low energy SO, S2 and P waves below 0.8 GeV, are somewhat displaced from what one 
gets by direct fits to experimental data, and fail to pass a number of consistency tests; notably, they do not 
fulfill well forward dispersion relations. All this is discussed in the present note, which is based on recent 
papers by the author and J. R. Pelaez, refs. 3, 4 and 5. 

2. The partial waves at low energy, s 1 / 2 < 1.42 GeV, from fits to data 



We will first consider wave-by-wave fits to experimental data for the SO, S2, P waves, as given in ref. 5. This 
improves the set of phase shifts, which was called "tentative solution" in ref. 3, and will provide us with a 
handy representation of what experiment gives for the nn scattering amplitudes. To fit the various phase 
shifts, ^ (/) (s), we use a parametrization for cot (^(s) which takes into account the analyticity properties of 
this quantity, as well as its zeros (associated with resonances) and poles (when the phase shift crosses mr, 
n =intcgcr). Explicit expressions will be given in each specific case. 



2.1. The SO, S2 and P partial waves at low energy, s 1 / 2 < 1 GeV 



We start with the P wave, for which we take the results obtained from a fit to the pion form factor, both in 
the spacelike and timelike regions, as given in ref. 6. We write 

cot (g) = ^(M 2 - s) (b + v^^ j , s V2 = 105GeV (21a) 

Zfc ^ -y/ S + y/ So — s J 

So is the point at which inelasticity begins to be nonegligible. We find 

So =1.069 ±0.011, Si = 0.13 ±0.05, M p = 773.6 ± 0.9, 
ai = (37.6 ± 1.1) x lO^M" 3 , 6j = (4.73 ± 0.26) x lO^M" 5 ; 



(2.1) 



this parametrization we take to be valid up to s 1 / 2 = 1 GeV. The fit to hit data is shown in Fig. 1. 

For the S2 wave we fit data where two like charge pions are produced: ^ although these pions are 
not all on their mass shell, at least there is no problem of interference among various isospin states. In the 
low energy region, we fix the Adler zero at z 2 = M n and fit only the low energy data, s 1 / 2 < 1.0 GeV; later 
on we will allow z 2 to vary. We have 

cot 5^ (s) = -^r~ — ~l-2 (go + B ^-^ZI l, Z2 = M„; sq 1 / 2 = 1.05 GeV, (2.2a) 



1 We use the convenient and self-explanatory notation of SO, S2, P, DO, D2, F, . . . , for the tytv partial waves. 
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FIGURE 1. The phase shifts of solution 1 from Protopopescu et al.,' 9 ' Hyams et al.' 10a ' and Es- 
tabrooks and Martin' 10c ' compared with the prediction with the parameters (2.1) (solid line below 
1 GeV). We emphasize that this is not a fit to these experimental data, but is obtained from the pion 
form factor. 161 The error here is like the thickness of the line. 

Above 1 GeV, the dotted line and error (PY) are as follows from the fit in ref. 5, Subsect. 2.1.2. 



B = -80.4 ±2.8, B x = -73.6 ± 12.6; 
(2) -- (-0.052 ± 0.012) M~\ b[ 2) = (-0.085 ± 0.011) M~ 3 



(2.2b) 



For SO, where the experimental situation is somewhat confused, we consider two different methods 
of data selection. In both, we fit those experimental data points in which on mass shell pions are involved: 
and K — > 2ir decays. In the first method, we also include some points at higher energy (0.81 GeV < 
5 1 / 2 < 0.97 GeV) where (most of) the various experimental numbers agree one with another within < 1.5 a. 
This we will call a global fit. Care is exercised to compose the errors in a realistic manner; the details may 
be found in ref. 5, Subsect. 2.2.2. In this case we fix the Adler zero at Zq = M w and find 



B =21.04, S x = 6.62, M a — 782 ± 24 MeV; x 2 / d -°-f- = 15.7/(19 - 3). 
a 0) =(0.230 ±0.010)M-\ 6 0) = (0.268 ± 0.011)M~ 3 ; 6§\m K )=4lJf ± 2.1 C 



(2.3a) 



this fit (shown in Fig. 2, PY curve) we take to be valid for s 1 / 2 < 0.95 GeV. The errors of the Bi are strongly 
correlated; uncorrelated errors are obtained if replacing the Bi by the parameters x, y with 

B = y-x; 51 = 6.62 - 2.590:; y = 21.04 ± 0.70, x = ± 2.6. (2.3b) 

Alternate possibilities are to fit only Kn and K — > 2tt data, or to add to this, individually, data 
from the various experimental analyses. The results can be found in Table 1 (see below, Subsect. 4.1). 
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FIGURE 2. The SO phase shifts and error band in the whole energy range as given by (2.3) below 
1 GeV, and from ref. 5 (Subsect. 2.2.4) above 1 GeV. The experimental points from Kn and de- 
cays are not shown. 

2.2. The SO, S2 and P waves at intermediate energy, 1 GeV < s 1 / 2 < 1.42 GeV 

We will not give here details on the fits to experimental data for the SO, S2 and P waves from ~ 1 GeV up 
to 1.42 GeV, which may be found in ref. 5. We will only comment that, while the S2 wave is reasonably well 
determined from experiment, the SO and P waves are not. Thus, for the first, depending on whether we use 
7T7T — ► 7T7T data, or data! 11 ! on 7T7t — > KK, the inelasticity 1 — 77^ varies by almost a factor three; we will 
discuss this again in connection with the scalar form factor of the pion in Sect. 9. Likewise, the values for 
<5i(s) given in the various experimental analyses disagree above 1.15 GeV. This means that relations, such 
as dispersion relations, in which the SO, P waves intervene, will have errors beyond the nominal errors above 
1 GeV. 

2.3. The D, F partial waves at low energy, s 1 / 2 < 1.42 GeV, from fits to data 

The D and F waves cannot be described with the same accuracy as the S, P waves. The data are scanty, 
and have large errors. To get reasonable fits we will impose the values of the scattering lengths that follow 
from the Froissart-Gribov representation. Note that this is not circular reasoning, and it only introduces a 
small correlation: the Froissart-Gribov representation for the DO, D2, F waves depends mostly on the SO, 
S2 and P waves, and very little on the DO, D2, F waves themselves. We do not discuss here the DO and F 
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FlGURE 3. Continuous line: The 1 = 2, D-wave phase shift, obtained by only fitting the ex- 
perimental data. Dashed line: with the parameters improved using dispersion relations. Dotted 
line: the fit, valid between s 1 ^ 2 = 0.625 GeV and 1.375 GeV, of Martin, Morgan and Shaw which 
ACGL and CGL, however, use from threshold to s 1 / 2 = 2 GeV. The experimental points are 
from Losty et al., Hoogland et al., solution A, and Cohen et al. (ref. 7) 



waves (that may be found in detail in ref. 5) as they do not present special features; only the D2 wave will 
be considered now. 

For isospin I = 2 we would only expect important inelasticity when the channel 7T7t — > pp opens up, 
so we will take the value sq = 1.45 2 GeV 2 ~ 4Af 2 for the energy at which elasticity is not negligible. A 

(2) 

pole term is necessary here to get an acceptable fit down to low energy, since we expect 5^ to change sign 
near threshold. The experimental measurements (Cohen et al.; Losty et al.; Hoogland et alJ 7 ') give negative 
and small values for the phase above some 500 MeV, while, from the Froissart-Gribov representation, it has 
been known for a long time' 12 ' that the scattering length must be positive. If we want a parametrization 
that applies down to threshold, we must incorporate this zero of the phase shift. The inflection seen in data 
around 1 GeV implies that we will have to go to third order in the conformal expansion. So we write 

cot4 2) - Sl/2r ~ ~ - ~ ^ 



>(-) - y {*o + B M s) + B M s?} w2 ;; 2) _ s (2.4a) 
with A a free parameter fixing the zero and 



w(s)= ^-^Z I g 1 / 2 = 1450 MeV . 

V s + V s - & 



Since the data we have on this wave are not accurate (cf. Fig. 3) we include in the fit the value of the 
scattering length that follows from the Froissart-Gribov representation, = (2.72 ± 0.36) x 10~ 4 M~ 5 . 
We get 

B a = (2.4 ±0.3) x 10 3 , B x = (6.8 ±0.8) x 10 3 , B 2 = (23.7 ± 3.8) x 10 3 , A = 196 ± 20 MeV . (2.4b) 

The fit (Fig. 3) returns reasonable numbers for the scattering length and for the effective range 
parameter, b^- 

4 2) = (2.5 ± 0.9) x 10~ 4 M^ 5 - b { 2 2) = (-2.7 ± 0.8) x 10~ 4 M/ 7 . (2.5) 
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3. The high energy (s 1 / 2 > 1.42 GeV) input 



In order to test the various dispersive sum rules we need an input for the imaginary part of the scattering 
amplitude at high energy (s 1 / 2 > 1.42 GeV) for which we here will take the standard Regge theory with the 
parameters of the various Regge poles as given in rcf. 4 (for the rho residue, we in fact take the improved 
version of ref. 5, Appendix B). We will not repeat here this Regge description but will note that, in the 
early 1970s, Pennington and Protopopescu (see e.g. rcf. 13) used crossing sum rules and then-existing low 
energy phase shifts data to shed doubt on the standard Regge picture for tttt scattering. 2 This was then 
accepted unqucstioningly by ACGL and CGL, who adopted a very distorted version of Regge behaviour 
with, for example, a Pomeron a third of what factorization (or experimental data on the total rnr cross 
section) implies, and unconventional slopes as well. This was used by ACGL, CGL in their analyses. As 
discussed in refs. 4, 5, however, standard Regge theory is perfectly consistent with crossing sum rules, as 
well as with experiment ' 14 ' if assumed to hold systematically above 1.42 GeV. 
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FIGURE 4. The -ktv cross sections. Experimental points from ref. 14. The stars at 1.38 and 
1.42 GeV (PY) are from the phase shift analysis of experimental data given in ref. 5. Continu- 
ous lines, from 1.42 GeV (PY): Regge formula, with parameters as in refs. 4, 5 (the three lines 
per fit cover the error in the theoretical values of the Regge residues). Dashed lines, above 2 
GeV: the cross sections following from ACGL;' 1 ' the grey band covers their error band. Below 
2 GeV, the dotted line corresponds to the ir + ir~ cross section from the Cern-Munich analysis; 
cf. Fig. 7 in the paper of Hyams et al.' 10 " 1 



2 It is appropriate to remark that such doubts could perhaps be maintained in 1974, when the tvtt phase shift were 
poorly known and, above all, it was also not clear which was the correct theory of strong interactions. They are of 
course difficult to believe once it became clear that the standard Regge picture is a feature of QCD. 
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In Fig. 4 we show the Regge values for the cross sections of the various irir scattering amplitudes, 
as given in refs. 4, 5, together with experimental data,! 14 ' a nd what follows from the imaginary parts used 
by ACGLW and CGLJ 2 ' It is not necessary to comment on this figure, which speaks for itself. 

We next remark that ACGL and CGL do not use Regge theory for 1.42 GeV < s 1 / 2 < 2 GeV, but 
take here the scattering amplitude as reconstructed from the phase shift analyses of the Cern-Munich group. 3 
Unfortunately, in this region inelasticty is large and the Cern-Munich experiments, which only measure the 
differential cross section for tttt — > tttt, are insufficient to reconstruct without ambiguity the full imaginary 
part. This is discussed in some detail in ref. 3 and, especially, in Appendix C of ref. 5 where it is shown that 
the Cern-Munich phases fail to pass a number of consistency tests. This is also seen very clearly in Fig. 4, 
where we plot the total cross section for 7r + 7r~ that follows from the analysis of Hyams et al.,[ 10a l which 
ACGL and CGL follow. Between s 1 / 2 <~ 1.5 and s 1 / 2 ~ 2 GeV, it is clearly incompatible with the other 
experimental data, I 14 ' as well as with Regge behaviour (as follows from factorization or fit to mr data). 

The situation is even worse for some specific waves. In particular, we will say a few words on the 
D2 wave. For D2, ACGL and CGL take (Appendix B.3 in ACGL) an empirical fit: 

4 2) (.s) = -0.003(s/4Af 2 ) (1 - 4M 2 / S ) 5/2 . (3.1) 

Eq. (3.1) was borrowed by ACGL from an old fit in the book of Martin, Morgan and Shaw, I 15 ' where 
only intermediate energy data were fitted. In fact, (3.1) fails at threshold (it gives a negative scattering 
length) and it does not fit well experimental data below 1.42 GeV, as shown in Fig. 3. Above 1 GeV, the D2 
wave in (3.1) grows (in modulus) quadratically with the energy, while from Regge theory it is easy to verify 
that all waves should go to a multiple of w and, in particular, D2 should go to zero; see Appendix C to ref. 5 
for details. It is true that this D2 wave is small but, given the accuracy claimed by CGL, it is certainly not 
negligible. 

It follows that the imaginary parts of the scattering amplitudes that ACGL, CGL use above 1.42 GeV 
are noticeably different from the physical ones. We will later see how much this affects their output. 



4. Improving the parametrizations with the help of dispersion relations 



4.1. Checking forward dispersion relations 



We expect that the scattering amplitudes that follow from the experimental phase shifts (and inelasticities) 
at low energy (s 1 / 2 < 1.42 GeV), together with the Regge expressions at high energy, will provide nn 
amplitudes that satisfy dispersion relations since they fit well the experimental data and are therefore good 
approximations to the physical scattering amplitudes. In the present Section we will check that this is the 
case, at low energies (s 1 / 2 5s 1 GeV), for three independent scattering amplitudes, which we will conveniently 
take the following i-symmetric or antisymmetric combinations, that form a complete set: 7r°7r° — ► 7r°7r°, 
7r°7r + — ► 7r°7r + , and the amplitude It — I, corresponding to isospin unity in the t channel. The reason for 
choosing these amplitudes is that the amplitudes for 7r°7r° and 7r°7r + depend only on two isospin states, and 
have positivity properties: their imaginary parts are sums of positive terms. Because of this, the errors are 
much reduced for them. 

We therefore consider here the following dispersion relations: for 7r°7r°, 

tj p p F MAf 2^_ *(*-4M 2 ) f 00 (2 S '-4M 2 )Im.F 00 ( 5 ') 

RcFooGs)-F 00 (4M 7r )_ P-P-^d- s ^ _ s)(s , _ 4 M 2 )( S ' + S - 4M 2 ) ' ^ 

In particular, for s = 2M 2 , which will be important when we later discuss the Adler zeros (Subsect. 4.2), we 
have 

a* <4.ib) 



For the SO wave ACGL take the re-elaboration of the Cern-Munich data by Au, Morgan and Pennington. 



[10d] 
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For 7r°7r+: 

At the point s = 2M 2 , this becomes 

F 0+ (4M 2 ) = F 0+ (2M 2 ) + D 0+ , D 0+ = f° ds ^Z°+(f) . (4 .2b) 

71" ./sm* s(s - 2M^)(s - 4M 2 ) 

Finally, for isospin unit exchange, which does not require subtractions, 

The imaginary parts at high energy, s 1 / 2 > 1.42 GeV, are taken from Regge theory as discussed in refs. 4, 5 
and in the previous Section. Depending on the experimental data we use to fit the SO wave we find the results 
reported in Table 1 for the fulfillment of the dispersion relations. The x 2 /d.o.f. is obtained by evaluating left 
and right hand sides of Eqs. (4.1a), (4.2a) and (4.3) at intervals of 25 MeV, from threshold to 0.925 GeV, 
adding the x 2 's for all points, and dividing by the number of points. We do not give, in Table 1, the x 2 /d.o.f. 
corresponding to 7r°7r + scattering, which is equal to 1.7 independently of the SO wave. 

We also give in Table 1 the central values of the quantity <5q ^ (0.8 2 GeV 2 ), that we will discuss later. 





Bo 


Si 


M a (MeV) 


dif.c- 1 ) 


V 2 

d.o.f. (7r " } 


4 0) (0.8 2 ) 


PY, Eq. (2.14) 


21.04 


6.62 


782 ± 24 


0.3 


3.5 


91.9° 


K decay only 


18.5 ± 1.7 


EE 


766 ± 95 


0.2 


1.8 


93.2° 


K decay data 
+ Grayer, B 


22.7 ± 1.6 


12.3 ±3.7 


858 ± 15 


1.0 


2.7 


84.0° 


K decay data 
+ Grayer, C 


16.8 ±0.85 


-0.34 ± 2.34 


787 ±9 


0.4 


1.0 


91.1° 


K decay data 
+ Grayer, E 


21.5 ±3.6 


12.5 ±7.6 


1084 ±110 


2.1 


0.5 


70.6° 


K decay data 
+ Kaminski 


27.5 ±3.0 


21.5 ± 7.4 


789 ± 18 


0.3 


5.0 


91.6° 


K decay data 
+ Grayer, A 


28.1 ± 1.1 


26.4 ±2.8 


866 ±6 


2.0 


7.9 


81.2° 


K decay data 
+ EM, s— channel 


29.8 ± 1.3 


25.1 ±3.3 


811 ±7 


1.0 


9.1 


88.3° 


K decay data 
+ EM, i-channel 


29.3 ± 1.4 


26.9 ±3.4 


829 ±6 


1.2 


10.1 


85.7° 


K decay data 
+ Protopopcscu VI 


27.0 ±1.7 


22.0 ±4.1 


855 ± 10 


1.2 


5.8 


82.9° 


K decay data 
+ Protopopescu XII 


25.5 ±1.7 


18.5 ±4.1 


866 ± 14 


1.2 


6.3 


82.2° 


K decay data 
+ Protopopescu 3 


27.1 ± 2.3 


23.8 ±5.0 


913 ±18 


1.8 


4.2 


76.7° 



PY, Eq. (2.3): our global fit, Eq. (2.3). We do not give errors here as they are strongly correlated. Grayer A, B, C, E: 
the corresponding solutions in the paper of Grayer et al.' 106 ' (note that solution D in this paper concerns only data 
at higher energies). EM: the solutions of Estabrooks and Martin. ' 10e ' Kaminski: the papers of Kaminski et al.' 10 ^' 
Protopopescu VI, XII and VIII: the corresponding solutions in ref. 9. 

Table 1 
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4.2. Improving the fits using forward dispersion relations 

In this Section we will improve the values of the parameters we have found with our fits to data requiring also 
fulfillment (within errors) of dispersion relations up to 0.925 GeV. This will provide us with a parametrization 
of the various waves with central values more compatible with analyticity and s — u crossing. This method 
is an alternative to that of the Roy equations to which it is superior in that we do not need as input the 
values of the scattering amplitude for \t\ up to 30M 2 , where the various Regge fits existing in the literature 
disagree strongly one with another (see ref. 5, Appendix B), and also in that, with dispersion relations, we 
can test all energies, whereas the Roy equations are only valid for s 1 / 2 < \/60 ~ 1.1 GeV (and in practice 
only applied up to 0.8 GeV). With the SO wave in Eq. (2.3), we find (units of M v ) 

SO; s 1/2 < 2m K : B = 17.4 ±0.5; Si =4.3 ±1.4; 

M a = 790 ± 21 MeV; z = 195 MeV [Fixed]; 
4 0) = 0.230 ± 0.015; b { ° } = 0.312 ± 0.014. 
S2; s 1/2 < 1.0 GeV : B = -80.8 ± 1.7; 5i = -77±5; z 2 = 147 MeV [Fixed]; 

4 2) = -0.0480 ± 0.0046; = -0.090 ± 0.006. 

S2; 1.0 < s 1/2 < 1.42 : 5 = -125±6; 5i = -119±14; e = 0.17 ±0.12. 

P; s 1/2 < 1.05 : B = 1.064 ± 0.11; B l = 0.170 ± 0.040; M p = 773.6 ± 0.9 MeV; 
ai = (38.7 ± 1.0) x 10~ 3 ; h = (4.55 ± 0.21) x 10~ 3 . 
DO; s 1/2 < 1.42 : B = 23.5 ±0.7; Si =24.8 ±1.0; e = 0.262 ± 0.030; 

4 0) = (18.4 ± 3.0) x 10~ 4 ; 4° } = (-8.6 ± 3.4) x 10~ 4 . 
D2; s 1/2 < 1.42 : B = (2.9 ± 0.2) x 10 3 ; B l = (6.3 ± 0.8) x 10 3 ; B 2 = (25.4 ± 3.6) x 10 3 ; 
A = 212 ± 19; 

4 2) = (2.4 ± 0.7) x 10" 4 ; bf> = (-2.5 ± 0.6) x 10" 4 . 

F;s 1/2 <1.42: B = (1.09 ± 0.03) x 10 5 ; B x = (1.41 ± 0.04) x 10 5 ; a 3 = (6.0 ± 0.8) x 10 -5 . 

(4.4) 

When using the fits to individual sets of data we get somewhat different results, for the SO wave parameters 
(the other waves are essentially like in (4.4), however): 



Improved fits: 


Bo 


Bi 


M a (MeV) 


z (MeV) 


x\h = 1) 
d.o.f. 


x 2 (oo) 

d.o.f. 


x 2 (o±) 

d.o.f. 


(4.1b) 


C(-8 2 ) 


PY, Eq. (4.4) 


17.4 ±0.5 


4.3 ±1.4 


790 ± 30 


198 ±21 


0.40 


0.66 


1.62 


1.6a 


91.3° 


K decay only 


16.4 ±0.9 


EE 


809 ± 53 


182 ± 34 


0.30 


0.29 


1.77 


1.5a 


91.3° 


K decay data 
+ Grayer, C 


16.2 ±0.7 


0.5 ±1.8 


788 ±9 


184 ± 39 


0.37 


0.32 


1.74 


1.5a 


91.0° 


K decay data 
+ Grayer, B 


20.7 ± 1.0 


11.6 ± 2.6 


861 ± 14 


233 ± 30 


0.37 


0.83 


1.60 


4.0 a 


85.1° 


K decay data 
+ Grayer, E 


20.2 ± 2.2 


8.4 ±5.2 


982 ± 95 


272 ± 50 


0.60 


0.09 


1.40 


6.0 a 


78.0° 


K decay data 
+ Kaminski 


20.8 ± 1.4 


13.6 ±3.7 


798 ± 17 


245 ± 39 


0.43 


1.08 


1.36 


4.5 a 


91.8° 



PY, Eq. (4.4): from the fit (2.3), improved with forward dispersion relations. Other conventions as in Table 1. Errors 
are given for the Adler zero, but we fix it at its central value when evaluating the errors for the other parameters. 

Table 2 
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5. The ACGL, CGL phases at the matching point, s 1 / 2 = 0.8 GeV 

A very important role is played, in the analyses of ACGL, CGL, by the input phases for the SO, S2 and P 
waves at the point, s 1 / 2 = 0.8 GeV, where they match the solutions of the Roy equations to the experimental 
amplitude at high energy. The errors these authors take for the S2, P waves are reasonable; that for the SO 
wave, which is a key quantity, is not, as we will now discuss. 

The quantity 5^ } {{0.8 GeV) 2 ) is in fact given in Eq. (6.3) of ACGL as 

4 0) ((0.8 GeV) 2 ) = 82.3 ±3.4°. (5.1) 

This value is also used in CGL. This error, ±3.4° , may be contrasted with the error estimates of ref. 5, 
Eqs. (2.13), which vary, for the data above 0.8 GeV, between 6° and 18°. 

Let us see if the ACGL number is reasonable. ACGL consider the difference Si — 5^ at 0.8 GeV, 
in the hope that some of the errors will cancel for this combination. They give a value for this quantity of 

4 O) ((0.8 GeV) 2 ) - tfi((0.8 GeV) 2 ) = 26.6 ± 2.8° . 

Now, the Cern-Munichl 10 ! experimental values for this which ACGL quote (there are others) are 

( 23.4 ±4.0° [Hyams ct al.] 
<5i((0.8 GeV) 2 ) - 4 0) ((°- 8 GeV) 2 ) = I 24.8 ± 3.8° [Estabrooks and Martin, s-channel] 

[ 30.3 ± 3.4° [Estabrooks and Martin, t-channel]. 




FIGURE 5. The SO phase shift corresponding to the improved fit in Eq. (4.4) (PY, thick con- 
tinuous line and error band), the unimproved solution of Eq. (2.3) (thin continuous line), and 
the improved solutions "K decay only" and "Grayer C" of Table 2 (difficult to see as they fall 
almost on top of PY). The solution CGL' 2 ' (lowest discontinuous line) is also shown. 
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FlGURE 6. Phase shift for the S2 wave, Eq. (4.4) (PY, thick continuous line, and error band); 
unimproved fit, Eq. (2.2) (PY, dotted line); and the phase CGL,' 2 ' and error band (dashed line). 

All the numbers here stem from the same experiment, and differ only on the method of analysis. 
Their spread is an indication of the systematic, or theoretical, error in the determination of ft — ^q '. To 
take this into account a correct procedure is to average the various determinations and enlarge the error, 
with what we will call "systematic error", so that all the numbers are covered by it: 

ft((0.8 GeV) 2 ) - 4° } ((0-8 GeV) 2 ) = 26.3° ± 3.8 (St.) ± 4.0 (Sys.). 

Using now ACGL's central value ft ((0.8 GeV) 2 ) = 108.9 ± 2° we would find 

^ 0) ((0.8 GeV) 2 ) = 82.7 ±7.8°; (5.2) 

we have added (St.) and (Sys.) errors linearly, as is mandatory when a single datum is involved. In fact, 
this is probably optimistic. First of all, the Hyams et al. value quoted above is only one of five solutions by 
the same experiment (depending on the corrections applied: cf. Grayer et alJ 10 ^). Secondly, if we included 
the data of Protopopescu et al./ 9 ! the error would increase to 10° . In any case, a realistic error should be 
several times what ACGL and CGL assume. 

But the best way to see that the error, and indeed, even the central value, given in (5.1) are not 
reasonable is to look at the experimental results reported in Tables 1, 2 for <5q ((0.8 GeV) 2 ), which speak 
for themselves. 

The character of a forced fit is also indicated in Fig. 5. CGL constrain the curve to pass through the 
point at (5.1). The corresponding solution is then dragged away from all solutions that are consistent with 
dispersion relations. In fact, it probably the distortion of the SO wave, together with the distortion caused 
by incorrect high energy behaviour, what also drives the S2 wave of CGL away from what one finds from 
experiment (Fig. 6). 

6. Dispersion relations and the CGL solution 

The inconsistency of the CGL solution with experimental results at high energy is particularly transparent 
if we consider the fulfillment of dispersion relations with the parameters of CGL for the SO, S2 and P waves 
al low energy, or with our parameters here. 
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FIGURE 7. Consistency of disper- 
sion relations for the 7T7t amplitudes of 
ref. 2 (CGL) and for our amplitudes, 
with the parameters in (4.4), denoted by 
PY. We plot the differences Aj, given in 
Eqs. (6.1), between the results of the cal- 
culation of the real parts directly with 
the various parametrizations, or from the 
dispersive formulas. Perfect consistency 
would occur if the continuous curves co- 
incided with the dotted lines. The error 
bands are also shown. The progressive 
deterioration of the CGL results as the 
energy increases is apparent here. 



This is depicted in Fig. 7, where we show the mismatch between the real part and the dispersive 
evaluations, that is to say, the differences Ai, 



it J AMI I s - s)(s' + s - 4M ~) 
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\ 
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s 
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and 

These quantities would vanish, = 0, if the dispersion relations were exactly satisfied. 

Wc include in the comparison of Fig. 7 the errors; in the case of CGL, these errors are as follow from 
the parametrizations given by these authors in ref. 2, for s 1 / 2 < 0.8 GeV. At higher energies they are taken 
from experiment via our parametrizations. By comparison, we show the same quantities for our best results 
in the present paper, that is to say, with amplitudes improved by use of dispersion relations, Eq. (4.4). In 
both cases we have taken the Regge parameters from Appendix B in ref. 5. 

7. Comparison of low energy parameters of CGL, DFGS, KLL, and here 



We here present, in Table 3, the low energy parameters as given in CGL, as well as in two other recent 
evaluations, that use the Roy equations, by Descotes et alJ 15 ', that we denote by DFGS, and by Kamihski, 
Lesniak and Loiseau' 15 ', denoted by KLL. This is compared with what one finds fitting experimental data, 
improved with dispersion relations, as reported in ref. 5 and repeated here, that we denote by PY. The 
mismatches between many of the parameters of CGL and PY are apparent in Table 3; not surprisingly, these 
mismatches affect mostly those parameters which are sensitive to high energy: b\ , a>^ , bSp ■ 





DFGS 


KLL 


CGL 


PY 


(0) 

% 


0.228 ± 0.032 


0.224 ±0.013 


0.220 ± 0.005 


0.230 ±0.015 


(2) 

% 


-0.0382 ± 0.0038 


-0.0343 ± 0.0036 


-0.0444 ± 0.0010 


-0.0480 ± 0.0046 


b (0) 




0.252 ±0.011 


0.280 ± 0.001 


0.312 ±0.014 


fc (2) 




-0.075 ±0.015 


-0.080 ±0.001 


-0.090 ± 0.006 


ffll 




39.6 ± 2.4 


37.9 ±0.5 


38.4 ±0.8 (xl0~ 3 ) 


bi 




2.83 ±0.67 


5.67 ±0.13 


4.75 ±0.16 (xl0~ 3 ) 


4 0) 






17.5 ±0.3 


18.70 ±0.41 (xl0~ 4 ) 








1.70 ±0.13 


2.78 ±0.37 (xl0~ 4 ) 








-3.55 ±0.14 


-4.16 ±0.30 (xl0~ 4 ) 


fc (2) 






-3.26 ±0.12 


-3.89 ±0.28 (xl0~ 4 ) 


as 






5.6 ±0.2 


6.3 ±0.4 (xl0~ 5 ) 



Units of M„. The numbers in the CGL column are as given by CGL in Table 2 and elsewhere in their text. 

In PY, the values for the D, F waves parameters are from the Froissart-Gribov representation. The rest are from the 

fits, improved with dispersion relations, except for oi and bi that have been taken as in ref. 5. 

Table 3 
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8. Comments on b\, on the scalar form factor of the pion 
and on chiral perturbation theory 



8.1. The b\ parameter, a sum rule and chiral perturbation theory 



The parameter b\ is one for which the different determinations in Table 3 vary more; for example, the CGL 
number is more than 4 a away from the result of PY, which follows from a very robust evaluation using the 
pion form factor in ref. 5. For this reason, we will look at it from two more points of view. First, we will 
write a sum rule for b\ that depends little on the high energy behaviour; and, secondly, we will look at b\ 
from the point of view of chiral perturbation theory. 

To get the sum rule we write a dispersion relation for the quantity 

0_ / F( f ' =1 )(a,0) \ 
ds \ s- 4M2 J ' 

which we evaluate at threshold. Taking into account that 

8_ / F( J ' =1 )( 8 ,0) \ = 3M^ 

ds { s -AMI J S =4MZ 2rr 15 

we obtain a fastly convergent relation for b\ : 



MM 



ds 



AMI 
1 



(s - 4M2)3 s 3 



ImF (7t - 0) (s,0) 



+ 



[s-AM^f s 3 



ImF (fl=1) (s,0) 



(s - 4M2)3 s 3 



ImF< /t=2 )(s,0) }. 



Most of the contribution to b\ comes from the SO and P waves at low energy, while all other contributions 
(in particular, the Regge contributions) are substantially smaller than 10~ 3 . Adding all pieces we find 

by = (4.99 ± 0.21) x 10~ 3 M~ 5 , (8.1) 

a value reasonably compatible with what we found in (4.4), with a completely different method, b\ = 
(4.55 ± 0.21) x 10~ 3 M~ 5 (because of correlations, the distance is actually ~ la). We can combine both 
and find a precise estimate, 

h = (4.75 ± 0.16) x 10" 3 M~ 5 . 

which is in fact the number reported in Table 3. 

We next turn to ch.p.t. To one loop we have! 17 ! 

M. 



1 288^ { ' h+h + ^ 



and /tt is the pion decay constant. If we replace the values of the k given in CGL we would find (in units 
with M n = l) 

h = (4.09 ± 0.81) x 10~ 3 [one loop; CGL constants h }. (8.2) 

In order to bring this into agreement with the CGL number, bi = (5.67 ± 0.13) x 10~ 3 , we would need 
corrections (two loop and higher) of <~ 40%. 

One may get an interesting relation for b\ by eliminating the constants l\ , li between b\ and a + , a o 
defined by 

a 0+ = |[4 0) -4 2) ], a o = |[4 0) +24 2) ]- 

We get 



h = | [3a 0+ - a 00 ] + + |) 



1 



288n 3 ttM w 



,3) 



- 13 - 



-F. J. YNDURAIN- 



Although this is valid only to one loop, the large higher order (two loop) correction due to the factor l// 4 
in (8.1) cancels almost completely. Moreover, the r.h.s. is dominated by ao+, which is known with great 
accuracy:' 5 ! one has ao+ = (10.61 ± 0.14) x 10 _4 M~ 5 . If we take the values of CGL for a o and a +, we 
obtain, from (8.3), 

h = (4.95 ±0.11) x 10~ 3 CGL. (8.4) 

This again is separated by several standard deviations from the value given in CGL. Contrarily to this, if we 
use the values in ref. 5 for the aoo, ao+, we find 

h = (4.44 ± 0.21) x 10~ 3 PY. (8.5) 

This number agrees with the value obtained with completely different methods (Table 3), thereby showing 
the consistency of the PY phase shifts, with themselves and with chiral perturbation theory. 

8.2. Chiral perturbation theory and the scalar form factor of the pion 

Finally, we compare the low energy 7T7t parameters with the results of chiral perturbation theory (ch.p.t.). 
We will consider the cases of CGL and our (PY) analyses. In both we take the value of the parameter 
F 3 = 2.9 ± 2.4 from the old calculation of Gasser and Leutwyler,! 17 ! although we will discuss it a bit more 
later. This parameter has very little influence in the results. We include in the comparison the quadratic 
scalar radius of the pion, (r| 7r ). For the CGL case we take it as given by the Donoghue, Gasser and Leutwyler 
calculation, I 18 l which CGL accept; for the PY case we take the determination obtained in ref. 19. This will 
be now discussed. 

Their results (from experimental analysis of the scalar pion form factor) for the pion scalar radius, 
(r§ „.), and of the ch.p.t. constant U, related to it by 

( r S,J = ^2f2 ~ if}' 

differs between the calculations of Donoghue et al.,! 18 ' who find U — 4.4 ±0.2, and those in ref. 19, where the 
value I4, — 5.4 ± 0.5 is obtained. Recently, Ananthanarayan et alJ 20 ' have repeated the calculation of ref. 18 
and conclude that it is confirmed (and thus that the result of ref. 19 is rejected by ~ 2 a). Unfortunately, 
the authors in ref. 20 still use the old parametrization of Hyams et al. (Eq. (12a) and Table 1 in ref 11a) 
for the K matrix for 7T7t and KK scattering to evaluate the scalar form factor of the pion, Fs(s). This is 
quite incompatible (by more than a factor three) with what one finds for the inelasticity, 1 — rj^ , by direct 
measurements of the tttt — ► KK cross section; I 11 ! seeFig. 8. Moreover, and unlike in the estimate of ref. 19, 
the authors in ref. 20 find a phase of the scalar form factor, S(s), clearly below 7r even for energies as high as 
s = 2 GcV 2 where from QCD estimates one expects it to be above it. In fact, by using methods like those 
for the vector form factor, 4 one finds that the scalar form factor of the pion behaves, at large s, as 

F s (s) ~ Cfl[ml{s)+fhl{s)]/s. 
Unlike for the vector case, however, the constant C can now not be calculated. From this it follows that 

Here A ~ 0.3 GeV is the QCD scale and d m = 12/ (33 — 2n f) is the anomalous dimension of the mass. Thus, 
and although part of the criticism of ref. 20 is correct in that the final state interaction theorem allows for 
a difference of tt in the phases of SO wave and form factor, we still consider that the result of ref. 19 is more 
reliable than that of refs. 18, 20: at least it is not incompatible with data on 7T7t — > KK scattering, and 

4 Note that he calculation is now less rigorous than what one has in the vector case. Details of the discussion on the 
scalar form factor of the pion, and the evaluation of I4, will be given in a forthcoming future article. 
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FlGURE 8. Fit to the 1 = 0, S- wave inelasticity and phase shift between 950 and 1400 MeV. Data 
from refs. 10a, lOf. The difference between the determinations of ?Jq from tttt — » iytt (PY from data) 
and from tt-r — > KK measurements' 11 ' (PY alternative) is apparent here. Note that the value of r/^ 
hat follows from the parametrization of Hyams et al.' lla ', that the authors in refs. 18, 20 use, lies at the 
lower part of the band denoted by "PY from data" here. 



agrees at high energy with QCD expectations. Moreover, it is consistent with what one finds fitting the 
experimental low energy parameters, as will be shown below. 

Another possibility would be to take the values for the chiral constant I4, from lattice calculations.' 21 ' 
The problem here is the existence of unknown, and potentially large theoretical and systematic errors. A 
last possibility is to obtain I4 by fitting the low energy parameters , b\, 6q > using for these last the ch.p.t 
results to one loop in ref. 17. We give here, for ease of reference, the various numbers one finds: 

4.4 ± 0.2 [refs. 18, 20] 
5.4 ±0.5 [ref. 19] 

4.0 ± 0.6 [lattice calculation, ref. 21] [s ' '> 

7.0±1.5exp.±0.11.o. [fitting 4 /} , b x , btf } ] 

(for the explanation of the errors "exp." and "l.o.", see below). 

The comparison of the low energy parameters and ch.p.t. for the CGL case is shown in Table 4; to 
get the one loop results there, we have taken the chiral lagrangian constants k as given by CGL themselves. 

It is revealing that almost all quantities require, to get agreement, rather large two loop corrections. 
This is suggestive of bias in the input, a bias that we have shown before to exist with the help of dispersion 
relations. 
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Quantity 


2 loop and other (CGL) 


1 loop ch.p.t. (CGL) 


l.o., ch.p.t. 


a (0) 




904 + 006 


0.157 


(2) 


— o D44-4- + n nm n 

VJ .\J L ± L ± L ± _L_ U.UU-LVJ 


—0 041 + 009 


—0.045 


b (o) 


n 980 + o oi o 


957 + 014 


0.179 




—0 080 ± 004 


—0 071 ± 007 


—0.089 


ffll 


(37.9 ±0.5) x 10~ 3 


(36.5 ± 0.8) x 10~ 3 


33.6 x 10~ 3 


bi 


(5.67 ±0.13) x 10~ 3 


(4.09 ±0.81) x 10~ 3 







(17.5 ±0.3) x 10~ 4 


(19.2 ±5.9) x 10~ 4 







(1.70 ±0.13) x 10~ 4 


(4.4 ± 1.6) x 10~ 4 







0.61 ± 0.04 fm 2 


0.57 ± 0.01 fm 2 






Table 4 



Comparison of evaluations of low energy parameters (according to CGL): leading order ("l.o.") and one loop (with 
the h parameters of CGL) and two loop chiral theory results and other results from CGL. 

We then turn to out to what one finds from experiment, i.e., using the low energy parameters from 
the 7T7T amplitude determined here and in rcf. 5, that we denote by (PY). We fix, from ref. 19, 

f 4 = 5.4±0.5, (8.7) 

and, from rcf. 17, 

f 3 = 2.9±2.5. (8.8) 

The D wave scattering lengths, and b\, may be written in terms of the pion decay constant / and 
the two chiral parameters l\, I2 only. For the first, 

(o) M" 1 
2 

(2) _ M- 1 n T 103 , (8.9) 



M- 1 
14407T 3 / 4 



" > - - - 7 — —{h + h — -jq- } ; 



the expression for bi was given in the previous Subsection. We have an ambiguity here: we may take the 
pion decay constant as equal to the physical constant, 

/ = U ^ 93.3 MeV; 

or we can identify it with the decay constant in the chiral limit, 



f M 2 - 1 

/=/o=/ n 1 "i6^% Z4 /- 



We will in fact take the average of both determinations, considering half the difference as an added error, 
called "h.o." , which will be a (rough) estimate of the corresponding higher order ambiguity. 

From the scattering length combination ao+ = §[c4°'' — a 2^]i which we determined very precisely 
(see Table 4 in ref. 5) we find 

h = 5.48 ± 0.06 (exp.) ± 0.61 (h.o.). (8.10a) 



From this and , 



l x = -0.76 ± 0.37 (exp.) ± 0.29 (h.o.). (8.10b) 
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In Eqs. (8.10), the first error is that due to the errors in the experimental and on l 3 , T 4 from (8.7), (8.8), 
and the second error is that due to the difference between using the two possibilities discussed for /, / = / 
or / = /„., in (8.9). Using this, we obtain the results of Table 5. 



Quantity 


Exp. value (PY) 


1 loop ch.p.t. 


l.o., ch.p.t. 


a (0) 


0.230 ± 0.015 


0.209 ±0.003 ±0.002 


0.157 


(2) 


-0.048 ± 0.005 


-0.042 ±0.001 ±0.001 


-0.045 


b (0) 

°0 


0.312 ± 0.014 


0.262 ±0.005 ±0.005 


0.179 


& (2) 


-0.090 ±0.006 


-0.073 ± 0.001 ± 0.003 


-0.089 


ai 


(38.4 ±0.8) x 10~ 3 


(37.6 ±0.6 ±0.6) x 10~ 3 


33.6 x 10~ 3 


bi 


(4.75 ±0.16) x 10~ 3 


(4.53 ± 0.21 ± 0.09) x 10~ 3 







(18.70 ±0.41) x 10~ 4 


input 







(2.78 ±0.37) x 10~ 4 


input 







0.77 ± 0.07 fm 2 


input 






Table 5 



Comparison of evaluations of low energy parameters from experiment (as determined in ref. 5), and from one loop 
chiral perturbation theory. Units of M„ 

For both and ftp -* already the one loop corrections (both for the CGL and PY evaluations) are 
quite large, so we expect also large two loop corrections, just by renormalization group arguments. In fact, 
a detailed estimate 1 22 1 for gives 

<5 2 ioo P ao 0) =0.017 ±0.002 
which brings the ch.p.t. value of, e.g., PY, of the SO scattering length to 

4 0) = 0.226 ± 0.003 ± 0.002 [incl. two loop], 

almost on top of the experimental value, 0.230 ± 0.015. For almost all other quantities, the two loop 
corrections necessary to bring the theoretical result in agreement with experiment are smaller than what 
CGL find; only the effective range parameters do not have overlapping error bars. 

We finish this article returning to the evaluation of Z3, I4. One can fix li, l 2 from the scattering 
lengths as in Eqs. (8.9), (8.10) and then fit I3, I4 from the all the other scattering lengths and effective 
range parameters. In this way one finds very poor x 2 /d.o.f.'s, which of course is just a reflection of the 
presence of higher order chiral corrections that we have not taken into account, 5 and the numbers 

F 4 = 7.1±0.7(exp.)±0.1(h.o.); k = 3.5 ± 10 (exp.) ± 2.9 (h.o.). (8.11) 

The value of I4 is almost identical to what we would obtain if fitting only &i, 6q 7 \ I4 — 7.0 ± 0.7 ± 0.1, that 
we had already reported in (8.7). This points clearly in the direction of what was found in ref. 19, and is 
well away from the results of refs. 18, 20. The value of I3 is compatible (within its huge error) with (8.8). 



5 The error denoted by "h.o." in our calculations only takes into account the estimated influence of the higher order 
corrections for the on Fi, F2. 



- 17 - 



-F. J. YNDURAIN- 

Appendix. A comment on a recent paper by Caprini, Colangelo, Gasser and 
Leutwyler 

In a recent paper Caprini, Colangelo, Gasser and Leutwyler, I 23 l to be denoted by CCGL, review the work of 
Pelaez and the present author in ref. 3 and conclude that, still, they consider the CGL solution consistent. 
We have already discussed this in connection with the dispersion relations, from which one can conclude that 
the claims in CCGL are unjustified. We will here examine only their contention that the PY Reggeistics 
cannot be correct because it does not satisfy certain sum rules. Of course this contention is meaningless 
since, as shown in the text (cf. Fig. 4) the PY cross sections are perfectly compatible with high energy 
(s 1 / 2 > 1.42 GeV) experimental data, while the ACGL ones are not. However, there is perhaps some point 
in spending a few lines to examine this contention of CCGL. 

CCGL evaluate the scattering length and effective range combinations ao+ = §[&2^ — a 2^L a oo = 
and the like combinations &o+i ^oo with the Froissart-Gribov projection (as in refs. 3, 5) or 
with the Wanders's sum rules. They define the differences Aao+, A6o+, . . . , between the results of the two 
evaluations and find, if using the Reggeistics of refs. 3, 4, 

Aa 0+ = -0.11 x 1(T 4 , Aa 00 = 0.034 x 10~ 4 , 

(A 1) 

A6 0+ = -0.13 x 10~ 4 , Afeoo = 0.090 x 10~ 4 , 

in units with = 1; cf. Sect. 9 in CCGL. From the fact that these numbers are not zero (as crossing 
symmetry would imply), CCGL conclude that the Reggeistics in refs. 3, 4 is not in "equilibrium" with the 
low energy scattering amplitude. 

However, errors should be taken into account, something that CCGL apparently forget. If we 
include only those errors following from the Froissart-Gribov representations in PY, and of these only those 
coming from D, F waves (since the S, and some of the P wave contributions cancel out) and from the Regge 
parameters, we should replace the relations (A.l) by 

Aa + = (-0.11 ± 0.09) x 10" 4 , Aa 00 = (0.03 ± 0.09) x 10~ 4 , 
A6 0+ = (-0.13 ±0.07) x 10~ 4 , A6 00 = (0.09 ±0.08) x 10~ 4 . 

That is to say, all the numbers are less than or around one standard deviation off zero, with one exception. 
This is A&o+, where —0.13 is to be compared to 0.07, in fact, 1.8 a. That among four quantities one 
is slightly off, by 1.8 a as A6 + is in our case, is statistically reasonable: you expect deviations of that 
order about 20% of the time and you get one in four. Not bad, particularly as this occurs for a quantity 
where the contribution of the slope of the It — 2 exchange is large and (as stated several times in ref. 3) one 
cannot take very seriously discrepancies where the derivative of the isospin 2 exchange amplitude contributes 
substantially. This agreement (A. 2) is the more remarkable because the low energy amplitude was obtained 
from fits to low energy tttt scattering: while the dominant Pomeron and P 1 Regge parameters come from a 
totally independent source, high energy irN and NN scattering, via factorization. 

Then, CCGL go on (still in their Sect. 9) to claim that, for the ACGL, CGL asymptotics, the sum 
rules above hold to a remarkable degree of accuracy because they find, in this case, 



Aa 0+ = -0.006 x 10~ 4 , Aa 00 = -0.009 x 10 ~ 4 , 
A6 + = -0.007 x 10~ 4 , Afeoo = -0.03 x 10~ 4 . 



They consider this better than (A. 1,2), and thus justifying their choice of Regge parameters, at least in an 
"effective" manner. 

The consistency of (A. 3) is not too surprising: ACGL fit their Regge parameters by requiring 
consistency of crossing sum rules. On the other hand, CCGL should not be so happy that the Eqs (A. 3) 
arc fulfilled. The relations (A. 1-3) are such that the S waves, and partially the P waves, cancel out. So, 
they depend very crucially on the D waves. Since the ACGL amplitudes between 1.42 and 2 GeV are very 
distorted (as discussed in Sect. 4 here and Appendix C in ref. 5), in particular for the D2 wave (Fig. 3), it 
follows that the Regge parameters they find must also be very distorted. So they are. Alternatively, and 
since the comparison with experiment of the Regge formulas of ACGL has revealed their inadequacy (Fig. 4), 
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it follows that the lower energy piece (s 1 / 2 < 2GeV) of these authors has to be irrealistic, which it certainly 
is. 
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